Subextensive Scaling in the Athermal, Quasistatic Limit of Amorphous Matter in 

Plastic Shear Flow 
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We present the results of numerical simulations of an atomistic system undergoing plastic shear 
flow in the athermal, quasistatic limit. The system is shown to undergo cascades of local rear- 
rangements, associated with quadrupolar energy fluctuations, which induce system-spanning events 
organized into lines of slip oriented along the Bravais axes of the simulation cell. A finite size scaling 
analysis reveals subextensive scaling of the energy drops and participation numbers, linear in the 
length of the simulation cell, in good agreement with the observed real-space structure of the plastic 
events. 
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The recent years have seen an important number of 
numerical and theoretical studies of plasticity in amor- 
phous materials. The microscopic picture of plastic de- 
formations which emerges from these studies, however, 
is still incomplete... at best, fragmented. Numerical evi- 
dence that plastic deformation involves highly heteroge- 
nous displacements of molecules led, early on, to the con- 
cept of "shear transformation zones" , which are expected 
to play, for amorphous solids, the role of defects in crys- 
tals . [ll S Most theoretical models of plasticity rely 
on this idea and, following Eshelby, 0| on the expecta- 
tion that elementary shear transformations are associated 
with quadrupolar energy fluctuations. Theoretical works 
indicate that the existence of quadrupolar elastic fields, 
and the consequent long-range interactions between shear 
transformation zones is an important mechanism that 
can induce strain localization and fracture in amorphous 
materials. 0, 0, Although quadrupolar energy fluctu- 
ations have been observed in a numerical model of dry 
foams, ^ they have never been seen in molecular sys- 
tems. 

This line of research should be contrasted with the 
phase space interpretation of plastic deformation recently 
proposed by Malandro and Lacks, on the basis of the 
inherent structure formalism, [iflj These authors study 
shear induced changes in the potential energy landscape, 
and the consequences of such changes on the macro- 
scopic mechanical behavior of glasses. In order to isolate 
these effects, Malandro and Lacks consider the quasi- 
static deformation of an amorphous material at zero- 
temperature, a protocol which has been used since early 
numerical studies as a means to bypass intrinsic limi- 
tations of molecular dynamics algorithms. 1\ For small 
deformations, the system follows shear induced changes 
of a local minimum (inherent structure) in the poten- 
tial energy landscape. Elementary catastrophic events 
occur when the local minimum in which the system re- 
sides annihilates during a shear-induced collision with a 



saddle point. The deformation of an amorphous mate- 
rial thus involves a series of reversible (elastic) branches 
intersected by plastic rearrangements (see figure QJ. 

The inherent structure formalism provides a precise 
definition of an elementary plastic rearrangement, but 
several questions arise about the spatial organization of 
these transitions: Are plastic events related to shear 
transformation zones and quadrupolar energy fluctua- 
tions? Do they involve spatially localized dynamical 
structures? If not, how do these structures scale with 
system size? Conflicting answers to these questions can 
be found in the literature. From measurements of partic- 
ipation ratio, Malandro and Lacks indicate that the el- 
ementary rearrangements they observe are localized. 
Durian and coworkers, for a model of foam (athermal 
by construction), observe a power-law distribution of 
energy drops at small strain rates, but with a system- 
size independent cut-off, indicating that no scaling be- 
havior is to be seen, unless at a very specific point in 
the jamming phase-diagram. jl3| A contradictory view- 
point is supported by the observations by Yamamoto and 
Onuki of increasing lengthscale suggesting the emergence 
of delocalized events, and critical behavior in the low- 
temperature, low strain-rate limit. 

In this work, we study spatial organization of elemen- 
tary transitions between inherent structures in quasi- 
static shear deformation, using the soft-sphere interac- 
tion potential used by Durian. This choice was ini- 
tially motivated by our intent to avoid delocalized struc- 
tures. With this model we observe: (i) quadrupolar 
energy fluctuations and cascades of these during single 
transitions between inherent structures, (ii) elongated, 
crack-like events that span the whole shear cell. Next 
(iii) we measure the distribution of energy drops and par- 
ticipation ratios and show that the size of typical events 
scales linearly with the length of the simulation cell. We 
thus show that scaling behavior is to be found in the 
quasi- static limit, and that it is associated with a cas- 
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FIG. 1: Potential energy as a function of strain during quasi- 
static shear of a 200x200 system. 



cade of spatially correlated quadrupolar energy fluctua- 
tions, reminiscent of theoretical considerations. 

We perform numerical experiments using simple shear, 
or so called Lees-Edwards, boundaries. Particles interact 
through the soft sphere-potential (121] : 







dtj) dij < 1 



> 1 



where, dij — 2{\\x'i 



I / {D, + Dj), with D,, the di- 



ameter of particle i. fc, being the only energy scale in 
the problem is set to unity. In order to prevent crystal- 
lization, a binary mixture is used with: Dl ~ = 
Dl"-^,Nl/Ns = ^ i where NUNs) is the num- 

5 

ber of large (small) particles. The athermal, quasistatic 
shear algorithm consists of two parts. First the simula- 
tion cell is deformed by a small amount with particle po- 
sitions fixed in reciprocal space (i.e. fixed relative to the 
Bravais axes of the simulation cell), producing an affine 
deformation in real space. Next the potential energy of 
the system is minimized with the shape of the simula- 
tion cell held fixed resulting in corrections to the affine 
deformation. Physically, the quasistatic algorithm corre- 
sponds to a material which is being sheared in a much 
shorter time than the thermally induced structural relax- 
ation time, but a much longer time than any microscopic 

times: Tmicro. < < ^ < < Tstruc. • 

The initial sample is prepared with a standard conju- 
gate gradient minimization applied to a random state. A 
fixed area simulation cell with a packing fraction of 1.0 is 
used for all systems. This density, well above the random 
close packing limit, was originally thought to preclude 
the emergence of non-localized structures. We use a 
strain step of size 10"'' for all simulations. Results will 
be presented for three ensembles of systems (with sizes: 
= 12.5x12.5, 25x25, and 50x50) and also one single 
200x200 system. 

Proceeding with a discussion of the single 200x200 sys- 
tem, figure n shows the potential energy of the relaxed 
configurations as a function of the applied shear strain 
for a small interval of strain from .15 to .17. The curve 
is composed of continuous segments, broken up by dis- 
continuous drops. Malandro and Lacks have demon- 
strated that each discontinuity arises from the destruc- 
tion of a potential energy minimum induced by the im- 
posed shear strain. In agreement with them, we find that 
the system is microscopically reversible upon changing 
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FIG. 2: Potential energy and sum of the squares of the forces 
as the system progresses through the minimization algorithm 
during the event circled in figure Q above. 



the sense of the strain during the continuous segments, 
but becomes irreversible across the discontinuities which 
constitute the fundamental plastic events. 

Next we look at the energy relaxation during a typi- 
cal plastic event, which is circled in figure ^ The energy 
and sum of the squares of the forces during the conjugate 
gradient descent for this single energy minimization are 
shown in figure El In these plots, the horizontal axis 
represents the amount of progress through the conju- 
gate gradient algorithm. We checked in smaller systems 
that minimization via integrating the equations of steep- 
est descent yields similar curves. In this latter case the 
horizontal axis can be directly interpreted as time, but 
steepest descent cannot be used in large systems due to 
its intrinsic inefficiency. In steepest descent dynamics, 
the time derivative of the energy is precisely the sum of 
the squares of the forces, and this relation holds reason- 
ably well for our conjugate gradient trajectories. Fig- 
ure El shows plateaus in the energy, which correspond to 
configurations where the forces are small and hence are 
very close to being mechanical equilibria. These config- 
urations are quasi- equilibria: deficient equilibria, each of 
which allows for an escape into a new quasi-equilibrium 
with lower energy. The system cascades through a se- 
ries of quasi-equilibria of decreasing energy until finally 
arresting in a true minimum. 

The observation of these transitions suggests that a 
typical plastic event, as circled on figured might be de- 
composed into elementary sub- events. Figure shows 
the resulting change in potential energy which occurs 
during the first force peak (in figureEl). The quadrupolar 
pattern is apparent. It is the first time such a field has 
been observed in an atomistic simulation. We emphasize 
that the energy dissipation field shown in figure Oli corre- 
sponds to a single elementary sub-event and contrast this 
with the work of Kabla and Debregeas Q who observe 
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FIG. 3: a) The change in potential energy during the first 
force pealc in figure |5| The color scale is linear with pure 
white (black) corresponding to a local decrease (increase) in 
energy per unit area equal to: 5 x 10~*. The orientation of 
this quadrupole is precisely what one would expect based on 
the direction of the principle axes of the applied shear strain, 
b) The local relative displacements (defined as the displace- 
ment of a particle with respect to the average displacement 
of its neighbors) that result from the entire cascade circled in 
figure Only displacement vectors which are larger than .1 
times the maximum are drawn. The black line is a guide to 
the eye oriented along the oblique Bravais axis. The particles 
inside the circled cluster are those which move during the first 
peak in the forces from figure |5| and produce the quadrupolar 
pattern shown in a). 



such a quadrupolar field in a mechanical film model only 
after averaging over many plastic events. We observe 
these quadrupoles generically during the onset of other 
typical events like the one circled in figure ^ however, 
after the onset the situation becomes more complex as a 
cascade is initiated in which the system proceeds through 
a series of such elementary sub-events. The spatial signa- 
ture of each elementary sub-event is, of course, noisy and, 
as the cascade proceeds, several elementary sub-events 
may occur concurrently and overlap in space. These ef- 
fects contribute to an increasingly complicated energy 
dissipation field in which it is often difficult to disentan- 
gle elementary quadrupolar patterns. Our observations, 
however, are consistent with viewing every cascade as a 
superposition of quadrupolar fields, each associated with 
an elementary sub-event. 

Turning to the spatial organization of the sub-events 
into a cascade, the local relative displacements of each 
particle are drawn in figure where the particles as- 
sociated with the onset are circled. In this picture, non- 
affine rearrangements cluster along the oblique Bravais 
axis of the simulation cell. We observe similar crack-like 
patterns generically in plastic events, aligned preferen- 
tially along the oblique or non-oblique (horizontal) Bra- 
vais axes. These patterns are reminiscent of the shear- 
bands which are predicted by several mesoscopic mod- 
els of plasticity, and are expected to emerge from the 
interactions of local shear transformations mediated by 
quadrupolar fluctuations in the elastic fleld 0, 0, IE H- 
We emphasize, however, that the patterns we observe 
here are certainly distinct from persistent shear bands: 



they are transient events which occur during a single, in- 
finitesimal strain step. The location and orientation of 
these patterns fluctuate as the system is sheared, and by 
no means can we identify the emergence of any station- 
ary heterogeneous deformation field, as observed in some 
molecular simulations [l^ . 

The existence of non-localized dynamical structures 
is consistent with the observation by Yamamoto and 
Onuki, [llj in molecular dynamics simulations of glass- 
forming systems, of an increasing correlation length in 
the limit where first temperature and then strain rate go 
to zero. They claim that their data is consistent with the 
existence of critical behavior in this limit, but were un- 
able to access the putative critical point due to inherent 
limitations in the molecular dynamics algorithm. Our 
algorithm is locked precisely at the T ^ 0,7 limit, 
which enables us to perform finite-size scaling analysis at 
this point. Before proceeding, we must emphasize that 
we study here a different molecular model which was not 
expected to display non-localized structures. 10 More- 
over, we now show that the crack-like patterns we ob- 
serve are responsible for the emergence of specific types 
of scaling which were not identified in previous numerical 
works. 

Results for the steady-state distribution of energy 
drops and participation number are given in figure 0] for 
different system sizes. We first examine the distribution 
of participation numbers N. We find a clear increase 
of the average with the linear system size, (N) ^ L, 
and a corresponding shift in the distribution. This scal- 
ing is consistent with the assumption that the dominant 
events are system spanning faults of length L with a typ- 
ical transverse length scale a which does not depend on 
system size. We next look at the distribution of energy 
drops: it is well described by a power law with an ex- 
ponential cutoff at large events. The power seems to 
be slightly smaller for our largest system, and ranges 
from .7 to .5, which is in rough accord with earlier re- 
sults [T^fl^. However, as was the case with the partici- 
pation number, we observe dramatic system size effects, 
with the average energy and the cut-off increasing lin- 
early with L. This scaling is consistent with the idea 
that the energy dissipated during a plastic event scales 
like {E) Eq{N) EqL, where Eq is the elementary 
energy released per quadrupolar fluctuation. 

Our observations differ from previous claims found in 
the literature. The distribution of participation num- 
bers was studied by Malandro and Lacks for a three- 
dimensional model of a glass. Q These authors concluded 
that the average participation number became indepen- 
dent of system size for large systems, but their data only 
disfavors an extensive {L^) scaling. From our observa- 
tions, we expect that, in three dimensions, plastic events 
are likely to organize into fault planes: in such a case, the 
scaling would become (TV) ~ , which is also consistent 
with Malandro and Lacks' data. i3| Here we come to an 
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FIG. 4: a),b) Distribution of participation number and energy 
drops for systems of lengths, 12. 5, 25, and 50. Insets: Scaling 
of the average participation number and energy drop with 
system length. 



important point: the non-extensivity of the participation 
number does not mean that structures are locahzed, as 
sub-extensive, system-spanning, structures may emerge. 

Our conclusion also differs from that of Tewari et. al. 
[l^ , on two-dimensional models of foams. We believe this 
discrepancy originates from a subtle consequence of their 
use of finite strain rate simulations. In quasistatic shear, 
the respective timescales of plastic events and shear are 
completely separated. However, the crack-like patterns 
we observe result from cascading sub-events, as informa- 
tion propagates through the system; the plateaus in fig- 
ure 121 correspond to times when information propagates 
with little dissipation. At finite strain rate, the spatial 
development of plastic events and the overall deformation 
of the material occur concurrently: the plateaus of fig- 
ure m become slightly tilted due to the overall energy 
increase induced by the finite strain-rate. The crite- 
rion Tewari et. al. use to separate individual plastic 
events stipulates that energy should decrease monoton- 
ically during a single event; it may misidentify quasi- 
equilibria for true equilibria, thus precluding the com- 
plete identification of elementary events. A finite shear 
rate may thus artificially "break" single plastic rearrange- 



ments into several spurious sub-events which have no sim- 
ple interpretation in the energy landscape. 

In conclusion, we have presented results on an atom- 
istic system sheared in the athermal quasistatic limit. We 
demonstrated the organization, during cascades, of ele- 
mentary quadrupolar plastic zones into lines of slip ori- 
ented along the Bravais axes of the cell. We proceeded 
to perform a finite size scaling analysis which revealed 
a linear system size dependence which indicates that the 
faultlike patterns of energy fluctuations play a major role 
in the emergence of scaling behavior. The overall picture 
which emerges from our simulations thus explains and 
clarifies various controversial claims found in the litera- 
ture, providing a first step toward a unified view of amor- 
phous plasticity based on both the energy landscape and 
real space pictures. 
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